Persistence Problem in Two-Dimensional Fluid Turbulence 
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We present a natural framework for studying the persistence problem in two-dimensional fluid 
turbulence by using the Okubo- Weiss parameter A to distinguish between vortical and extensional 
regions. We then use a direct numerical simulation (DNS) of the two-dimensional, incompressible 
Navier-Stokes equation with Ekman friction to study probability distribution functions (PDFs) of 
the persistence times of vortical and extensional regions by employing both Eulerian and Lagrangian 
measurements. We find that, in the Eulerian case, the persistence-time PDFs have exponential tails; 
by contrast, this PDF for Lagrangian particles, in vortical regions, has a power- law tail with an 
exponent = 2.9 ± 0.2. 
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The persistence problem, which is of central im- 
portance in nonequilibrium statistical mechanics pQ, 
is defined as follows: For a fluctuating field </>, the 
persistence-time probability distribution function (PDF) 
P^(r) yields the probability that the sign of qb at a point 
in space does not change up to a time r. Theoretical, 
experimental, and numerical studies of a variety of sys- 
tems, ranging from reaction-diffusion systems to granular 
media, have shown that P^(r) ~ r~ e as r — » oo, where 
is the persistence exponent This nontrivial exponent 
cannot be obtained from dimensional arguments; it can 
be calculated analytically only for a few models [2j; for 
most models it has to be obtained numerically. We pro- 
pose a natural way of defining the persistence problem 
in two-dimensional turbulence. We then show how to 
obtain the persistence exponent for this case. 

Turbulent flows in two-dimensional fluid films dis- 
play vortical points and strain-dominated or extensional 
points. We show how to examine the persistence of such 
points in time by direct numerical simulations (DNSs) 
of the forced, two-dimensional, incompressible Navier- 
Stokes equation. Our study has been designed with thin- 
fluid- film experiments in mind [3j [4] so we account for an 
air-drag-induced Ekman friction and we drive the fluid 
by using a Kolmogorov forcing. We demonstrate that the 
Okubo- Weiss parameter [5[ [6] A, whose sign at a given 
point determines whether the flow there is vortical or 
extensional, provides us with a natural way of studying 
such persistence. 

It is important to distinguish the following types of 
persistence times: (A) In the Eulerian framework we con- 
sider a point (x, y) and, by following the time evolution of 
A, determine the time r for which the flow at this point 
remains vortical (extensional) if the flow at this point 
became vortical (extensional) at some earlier time; (B) 



in the Lagrangian framework we consider how long a La- 
grangian particle resides in a vortical (extensional) region 
if this particle entered that vortical (extensional) region 
at an earlier time. For all these cases we obtain PDFs of 
the persistence or residence times that we denote gener- 
ically by r. We find, in the Eulerian framework, that 
the PDFs of r show exponential tails in both vortical 
and extensional regions. In the Lagrangian framework 
the PDF of the residence time of the particle in exten- 
sional regions also shows an exponential tail; the anal- 
ogous PDF for vortical regions shows a power-law tail. 
The persistence exponent that characterizes this power 
law is independent of parameters such as the Reynolds 
number, the characteristic scale of the forcing, and the 
coefficient of Ekman friction, at least at the level of our 
numerical studies. 

We perform a direct numerical simulation (DNS) of the 
incompressible, two-dimensional, Navier-Stokes equation 

d t uj - J(ip, uj) = vV 2 uj + fuj - fiuj, (1) 

with periodic boundary conditions, by using a pseu- 
dospectral method [7] with TV 2 collocation points and the 
2/3 dealiasing rule. Here ip is the stream function, uj the 
vorticity, J(ip,uj) = (d x tlj)(dyUj) — (d x uj)(dyip)^ v the kine- 
matic viscosity, and \i the coefficient of Ekman friction. 
At the point (x,y) the velocity u = (—d y ip^d x ip) and 
the vorticity uj = V 2 ^; the external deterministic force 
f UJ (x 1 y) = — Fo&inj cos(/ci n jx), where Fo is an amplitude 
and fcinj the energy-injection scale in Fourier space. The 
injected energy displays an inverse cascade to small k. 
Ekman friction removes energy from all Fourier modes; in 
particular, it removes energy from small-/c Fourier modes 
in such a way that the system reaches a nonequilibrium 
statistically steady state. We evolve Eq. [I] in time by a 
second-order, exponential Runge-Kutta method [8]. In 
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all our simulations we wait for a time T tran [see the cap- 
tion of Table Q] to allow transients to die out so that 
our system reaches a statistically steady state. Figure ([!]) 
shows a typical pseudocolor plot of ip in such a state. 

To calculate Lagrangian quantities we track N p parti- 
cles. The evolution equation for a Lagrangian particle is 
^x^(t) = u(x£,£), where x^(t) is the position of a La- 
grangian particle at time t; u(xl,£), the velocity at the 
Lagrangian particle position, is evaluated from the Eule- 
rian velocity field u(x, t) by using a bilinear-interpolation 
scheme [9]. The evolution equation of the particles is 
solved by a second-order, Runge-Kutta method [9 . Ini- 
tially all the particles are seeded randomly into the flow. 
A typical Lagrangian-particle track superimposed on a 
representative pseudocolor plot of ip is shown in Fig. ([!]). 
The list of parameters used in our simulations is given in 
Table @. 

From the velocity-gradient tensor A, with components 
Aij = diUj, we obtain the Okubo- Weiss parameter A, 
the discriminant of the characteristic equation for A. If 
A is positive (negative) then the flow is vortical (exten- 
sional) [5]. In an incompressible flow in two dimensions 
A = det^4; and the PDF of A has been shown [4] to 
be asymmetrical about A = (vortical regions are more 
likely to occur than strain-dominated ones). For the Eu- 
lerian case, we monitor the time evolution of A at N p 
randomly chosen points that are fixed in the simulation 
domain. In our Lagrangian study, we begin with the val- 
ues of A on the spatial grid that we use for our Eulerian 
DNS; bilinear interpolation then yields A at the positions 
of Lagrangian particles, which can be at points that do 
not lie on this grid; we can thus monitor the evolution of 
A along Lagrangian-particle trajectories. 

We denote the presistence-time PDFs by P and the as- 
sociated cumulative PDFs by Q; the subscripts E and L 
on these PDFs signify Eulerian and Lagrangian frames, 
respectively; and the superscripts + or — distinguish 
PDFs from vortical points from those from extensional 
ones. To find out the persistence-time PDF P# (t) [resp., 
Pg(r)] we analyse the time-series of A obtained from 
each of the N p Eulerian points and construct the PDF of 
the time-intervals r over which A remains positive (resp., 
negative). The same method applied to the time series 
of A, obtained from each of the N p Lagrangian particles, 
yields P+(r) [resp., P-(t)]. 

We use the rank-order method [10] to calculate cu- 
mulative PDFs because they are free from binning er- 
rors. We show representative plots of the cumulative 
PDFs Qe (red crosses), (black open circles), and 
(magenta full circles) in Fig. ; the dashed lines in- 
dicate exponential fits to these cumulative PDFs. From 
these and similar fits we conclude that the PDFs p£ , 
P^, and Pj7 have exponentially decaying tails, for all 
the runs, from which we can extract characteristic time 
scales. In particular, from the Eulerian PDFs we ob- 
tain the times and T^, for vortical and strain- 



dominated regions, respectively, which are defined as fol- 
lows: P+(r) - exp(-r/T+) and P~ (r) - exp(-r/T E ") 
as r — » oo [see Table (□)]. The Lagrangian PDF P^ also 
has an exponentially decaying tails from which we obtain 
the characteristic time [see Table 

The persistence-time PDF P^ and the associated cu- 
mulative PDF Q \ of a Lagrangian particle in a vortical 
region is very different from those discussed above: the 
tails of P^ and Q \ have power-law, and not exponential, 
forms; we show this in Fig. ^p) via a representative plot 
Q\ for the run R4. Thus, as in nonequilibrium statisti- 
cal mechanics [1 , we can define the persistence exponent 
via ~ r~ e . Our run R4 has the largest value of 
Re\ amongst the runs Rl — 4 [Table (|l|)] and, therefore, 
is best suited for estimating from plots such as the one 
in Fig. ([2J)). We obtain the exponent (—0 + 1) by fitting 
a power-law to the tail of the cumulative PDF Q\. To 
find the best estimate for (—0 + 1) we evaluate the local 
slope x = dlog 10 Q^(r)/dlog 10 (r) in the region shown 
in the inset of Fig. ^p). Our estimate for (—0 + 1) is the 
mean value of x over the region indicated in the inset; the 
standard deviation of x yields the error; finally we obtain 
= 2.9 + 0.2. Our other runs Rl — 3 yield smaller scal- 
ing ranges than the one in Fig. (J^Jd) and, therefore, yield 
values for with larger error bars; but these values are 
consistent with our estimate for from run R4. We also 
find that the persistence exponent does not depend on 
the parameters /i, Po, and fci n j (within error-bars) in the 
range of parameters accessible in our simulations. Based 
on this evidence we conjecture that is a new univer- 
sal exponent that characterizes two-dimensional, Navier- 
Stokes turbulence. A conclusive proof of this conjecture 
must await future studies. 

The exponent > 1, so we can obtain the average life- 
time T^ean °f a particle in a vortical region from P^; 
Another estimate of the lifetime of vortices in the La- 
grangian frame is given by the time scale T c f, which is 
the cutoff scale of the power-law decay of tail of P^, 
[Table ft}]. 

In most persistence problems in nonequilibrium sta- 
tistical mechanics the power-law tail for the presistence- 
time PDF appears in the following way: Typically we 
consider the persistence-time PDF of <\> that comes from 
a Gaussian, but non- stationary, process; a transforma- 
tion to logarithmic time s transforms such a process to a 
Gaussian stationary process (GSP) X(s); if, in addition, 
the GSP is also Markovian, i.e., the autocorrelation func- 
tion of the GSP, f(t) = (X(s)X(s+t)) s , is an exponential 
function of t, then the presistence-time PDF of the GSP 
can be shown to have an exponential tail pQ. If we now 
transform back from logarithmic to linear time, the expo- 
nential tail is transformed to a power-law tail. Further- 
more, for a GSP for which the autocorrelation function 
f{t) decays faster than 1/t for large t, the persistence- 
time PDF P(t) ~ exp(— 7^), where 7 is a constant [1 . 

In the two-dimensional, fluid-turbulence problem that 
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we study here, we have checked numerically that A is a 
stationary process; hence we do not need to transform 
to logarithmic time. We have calculated two types of 
autocorrelation functions for A (we denote these generi- 
cally by C\(t) in Fig. (J^)): (a) For the first we evaluate 
(A(O)A(t)) over the track of a Lagrangian particle; this is 
shown by blue open circles in Fig. (J^); (b) for the second 
we evaluate (A(O)A(t)) at a given point on our Eulerian 
grid; this is shown by red full circles in Fig. (J^); here (•) 
denotes averages over different origins of time and also 
over N p different Lagrangian particles, for case (a), or 
over TVp different Eulerian positions, for case (b). As we 
show in the inset of Fig. (J^) for both these cases, C\(t) 
is approximated well by the function exp[— (t/T\) 2 ], over 
the range 10 -4 < (t/Trj) < 10 _1 ; this decay is clearly 
faster than 1/t for large t. Here Ta, the characteristic de- 
cay time, is slightly larger in the Lagrangian case than in 
the Eulerian one; however, in both these cases Ta ~ T v . 
Plots of Ca from our other runs Rl — 3 are similar to the 
one in Fig. (J^) from run R4. Note that, in the problem 
we study, the persistence-time PDF is not constrained to 
have an exponential tail because A, although stationary, 
is not a Gaussian process. 

We have presented a natural framework for studying 
the persistence problem in two-dimensional fluid turbu- 
lence. The most important result of our study is that the 
persistence-time PDFs for vortical points show qualita- 
tively different behaviors in the Eulerian and Lagrangian 
cases: In vortical regions, for the Eulerian case, this PDF 
displays an exponential decay; in contrast, for the La- 
grangian case it shows a power- law tail. Qualitatively 
such nontrivial behavior appears because a passive par- 
ticle can be trapped for quite some time in a vortical 
region [11]) Furthermore, we provide a way of measur- 
ing the lifetime of a vortex precisely. In the Eulerian 
frame the characteristic lifetime of a vortex is the time 
scale that follows from the exponential form P^(r) ~ 
exp(— r/Tg"). However, in the Lagrangian frame the 
persistence-time PDF for vortical points shows a power- 
law decay with a persistence-time exponent 6 = 2.9 ±0.2. 
Hence there is no single time scale which describes the 
time spent by passive particles in vortical points; but, as 
we have mentioned above, an average residence time can 
be defined because > 1. 

The PDF of residence times of passive tracers in vor- 
tical regions P^ is of great fundamental and engineering 
importance. Earlier studies [I2j [13] have attempted to 
measure this PDF; however, their methods of obtaining it 
are not as precise as the one we present here: the Okubo- 
Weiss parameter A, which we employ, helps us to dis- 
tinguish clearly between vortical and extensional regions 
in the two-dimensional flows we consider. The natural 
way of generalizing our study to its three-dimensional 
counterpart is to replace the Okubo- Weiss parameter by 
QR plots [14 and then to study the PDFs of residence 
times of Lagrangian particles in each quadrant of the QR 
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FIG. 1: (Color online) (a) A pseudocolor plot of the stream 
function ip, at a representative time in the statistically steady 
state, with a representative Lagrangian particle track (blue 
squares) superimposed on it from our run R2. The symbol 
o indicates the beginning of the trajectory and the x sign 
marks its end. For an animated version see the movie file at 
http : //www . youtube . com/watch? v=ANDUNDHYgzk ; (b) log- log 
(base 10) plot of the energy spectrum for our run R4 (line with 
dots); the black line with a slope —3.6 is shown for reference. 



plot; here Q = -\tr(A 2 ) and R = — |tr(^4 3 ) where A 
is the velocity gradient matrix. It remains to be seen 
whether the tails of such PDFs have power-law tails in 
three-dimensional turbulence; such a study lies beyond 
the scope of this paper. From the point of view of the 
general theory of persistence problems, the time-series of 
A in the Lagrangian frame provides a particularly inter- 
esting example. The persistence-time PDF has a power- 
law behavior for positive A but an exponential tail for 
negative A. A similar, but less dramatic, example of 
such asymmetry has been observed for the case of grow- 
ing interfaces whose height h obeys the Kardar-Parisi- 
Zhang (KPZ) equation [15]. For this KPZ case, the 
persistence-time PDF shows power-law tails, but with 
different persistence-time exponents for positive and neg- 
ative h. We hope our study will stimulate new experi- 
mental investigations of persistence-time PDFs in two- 
dimensional fluid turbulence and also studies of such 
PDFs for other non-Gaussian, stationary processes. 
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Run N v /jl F fejnj Id A Re\ T eddy T v T£ T+ ean r in j T cf 

Rl 512 0.016 0.1 45 10 2.3 x 10" 2 0.2 59 0.1 3.4 x 10" 2 0.3 ± 0.04 0.12 ± 0.02 0.21 ± 0.05 2.8 x 10" 2 0.3 0.9 

R2 512 0.016 0.45 45 10 2.1 x 10" 2 0.1 27 0.1 2.7 x 10" 2 0.4 ± 0.05 0.17 ± 0.02 0.24 ± 0.03 4.2 x 10" 2 0.2 0.8 

R3 1024 10" 5 0.01 0.005 10 4.3 x 10" 3 0.1 827 11 1.9 19 ± 3 10 ± 2 13 ± 2 1.8 19.9 76 

R4 1024 10" 5 0.01 0.005 4 5.4 x 10" 3 0.2 1319 7 2.9 31 ± 5 15 ± 2 25 ± 4 2.5 30.2 81 



TABLE I: Parameters for our runs Rl-4: N is the number of grid points along each direction, N p = 1000 is the number of 
Lagrangian particles and Eulerian positions (at which we monitor A), v the kinematic viscosity, fi the Ekman friction, Fo the 
forcing amplitude, fa n j the forcing wavenumber, Id = (u 3 /e) 1 ^ 4 the dissipation scale, A = ^vE/s the Taylor microscale, Re\ = 
w rms A/^ the Taylor-microscale Reynolds number, T e dd y = [tt ^2 k (E(k)/k)/(2u^ ms )] the eddy-turn-over time, and T v = \fvjs 
the Kolmogorov time scale. The time scales TjjT, T£~ and T^" are obtained from exponential fits to the tails of the cumulative 
PDFs , Ql and , respectively, as shown in Fig ([^l) . T+ ean is the average time spent by a Lagrangian particle in a vortical 
region, Ti n j = (/ 2 n j/^i n j) 1//3 is the energy-injection time scale, where E\^ =< f u • u >, (f^ = V x f u ), is the energy- injection 
rate and /i n j = 27r//ci n j is the energy- injection length scale. T c f is the large-time cutoff of the scaling range of as shown in 
Fig.|2]3. We do not use data from the initial period of duration T tra n = 100T e dd y ; this removes the effects of transients. We use 
a square simulation domain with side L = 2tt and grid spacing S x = L/N. 
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FIG. 2: (Color online) (a) Representative semilog plots (base 10) of the cumulative persistence-time PDFs (red crosses), 
(black open circles), and (magenta full circles); the subscripts E and L signify Eulerian and Lagrangian frames, respectively; 
and the superscripts + or — distinguish PDFs from vortical points from those from extensional ones. The dashed lines indicate 
exponential fits to these cumulative PDFs. These plots use data from our run R4. (b) A representative log-log plot (base 10) 
of the cumulative PDF (r)(«) versus r for our run R4; the full red line, with a slope equal to —2, is drawn for reference. 
The vertical arrows indicate the time scales (from left to right) T v , Ti n j, and T c f, respectively. The inset shows the local slope 
X = dlog 10 (t) / d\og 10 (r) versus r; the horizontal line is drawn at the mean (x) — —1.93 and the black dashed lines, drawn 
at (x) =t cr x , where a x c± 0.2 is the standard deviation of %. (c) Plots of the autocorrelation function of the Okubo- Weiss 
parameter Ca^/T^) (see text) versus t/T^ in Eulerian (blue open circles) and Lagrangian (red full circles) frames for our run 
Rl. In the inset we compare the data points for these plots with their fits (full lines) to the form Gexp[— (£/Ta) 2 ]. 
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